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ABSTRACT 

We investigate the emergence of hyperons in black-hole-forming failed super- 
novae, which are caused by the dynamical collapse of nonrotating massive stars. 
We perform neutrino-radiation hydrodynamical simulations in general relativity 
adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive 
cases are examined for the potential of S hyperons. Since hyperons soften the 
EOS, they shorten the time interval from the bounce to black hole formation, 
which corresponds to the duration of neutrino emission. This effect is larger 
for the attractive case than the repulsive case because S hyperons appear more 
easily. In addition, we investigate the impacts of pions to find that they also 
promotes the recollapse towards the black hole formation. 
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Introduction 



The gravitational collapse of massive stars is the fate of the stellar evolution, which leads 
to the variety of explosive phenomena, the formation of compact objects and the neutrino 
burst. The type lb, Ic and II supernovae are thought to be driven by the core collapse due 
to photodisintegration reaction of nucleus. The collapse is bounced by the nuclear repulsion 
and the shock wave is launched invoking explosion. They lead to the neutron stars. In 
somewhat low mass cases, the core collapse may be triggered by electron captures before 
Ne ignition, which is called electron-capture supernovae. However, the mass range of these 
supernova prog enitors is still un c ertain . If the stellar mass is lower, a white dwarf will be 
formed finally. iPoelarends et al.l (120081 ) investigated stellar evolution sequences with initial 
masses between 6.5 and ISMq and the solar metallicity using three diffe rent codes and found 
that the lower mass limit for supernovae is 9-12M0. Observationally, ISmartt et al.l (|2009[ ) 
found that it converges to 8 ± IMq from the direct detections of progenitors. 

Massive stars beyond a certain mass limit lead to different fates from ordinary super- 
novae though the upper mass limit is more uncertain. ISmartt et al.l (120091 ) suggested that 
the majority of massive stars above ~2OM0 may collapse quietly to black holes and that 
the explosions remain undetected. They investigated 92 samples which complete the core 
collapse supernovae occurred in a fixed 10.5-year period within a distance of 28 Mpc and 
concluded that there in no evidence o f massive progenitors . Apart from the direct detections, 
according to the lightcurve models, iNomoto et al.l (120061 ) proposed that the massive stars 
beyond 25M0 form black holes and their fate splits into two branches, namely, a hypernova 
branch and a faint-supernova branch. While the constituents of the hypernova branch, which 
have possible relevance with gamma ray bursts, are suggested to be strongly rotating, those 
of the faint-supernova branch are nonrotating or weakly rotating. 

Having observational suggestions, a quiet death through the core collapse occurs for a 
certain fraction of the massive stars involving collapse. They are thought to become black 
holes eventually. Their fate is theoretically proposed to split further into two categories: 
prompt and delayed black hole formations. While the delayed black hole formation occurs 
after the faint (perhaps even undetected) supernov a explo s ion, t he prompt black hole for- 
mation occurs when the supernova explosion fails. iFryerl ( 119991 ) showed that nonrotating 
stars with the initial r nass of >40Mf7i b e come failed supernovae according to his numeri- 
cal simulations. Lately, lO'Connor fc OttI (120 111 ) performed the core collapse simulations for 
many progenitors with the spherically symmetric models involving simplified neutrino trans- 
fer. They found that the threshold between the prompt and delayed black hole formations 
depends severely on the progenitor models, or evolutionary calculations. 



For the investigation of the black hole formation, implementation of properties of dense 
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matter at supranuclear density is mandatory. In particular, meson condensation, hyperon 
appearance and quark deconfinement are expected, affecting the equation of state (EOS) 
markedly. For the delayed black hole formations, a proto-neutron star is formed once and 
recoUapse to a black hole >10 sec after the bounce. iKeil &: Jankal (Il995l ) investigated evolu- 
tionary calculations of the proto-neutron star cooling taking into account hyperons. Their 
result is adopted as the initial condition of the dyn amical simulation for th e final 100 ms 
before the delayed black hole formation performed by lBaumgarte et al.l ( 1l996l ). I ncidentally, 



the ev olutions of proto-neutron s tars are studied a lso with kaon condensation ( jPons et al. 



2001al ) and quark deconfinement (jPons et al.ll2001bl ). Those results clarified that the exotic 
phases appear at late stage of the evolution of proto-neutron stars for ~20 sec. 

The prompt black hole formation is a promising branch to explore the exotic phase, since 
the exotic phase appears dramatically right after th e core bounce. The prompt black hole for- 



mation occ urs within ~1 sec after the bounce (e.g.. lLiebendorfer et al.ll2004 



Sumiyoshi et al. 



20061 . |2007| ) depending on the EOS. While no electromagnetic signal other than the disap- 
pearance o f progenitors is expecte d for the failed supernova involving the prompt black hole 
formation ( iKochanek et al.l 120081 ) , it is as bright in neutrino emissions as ordinary core- 
collapse supern ovae. Recently, nume rical simulations of failed supernovae utilizing the EOS 
with hyperons fishizuka et al.l 120081) we r e perfo rmed to evaluate their neutrino signal by 
Sumiyoshi et al. ( 20091 ). Nakazato et al. ( 2010b) studied f urther the neutrino signal utiliz- 
ing the EOS with quarks and pions ( iNakazato et al.ll2008l ). Since these exotic constituents 
soften the EOS and reduce the maximum mass of neutron stars, the time interval between 
the bounce and black hole formation, which corresponds to the duration of neutrino emis- 
sion, gets shor ter. Therefore the neu trino signal can be used to diagnose the emergence of 
exotic matter ( iNakazato et al.ll2010al ). 



We carry forward here the investigation of the failed supernovae and prompt black hole 
formations with hyperons by paying attention to the ambiguity of hyperon interaction. The 
hyperon interactions with nucleons are relatively well known comparing with other exotic 
constituents (kaons, pions and quarks). In particular, A-A^ interaction is well investigated 
through the single-particle energies of A hypernuclei. Unfortunately, however, T,-N inter- 
action, which affects the components of dense matter and t he stiffness of EOS, h as a large 
uncertainity even at present. Within this S-A^ ambiguity, llshizuka et al.l (|2008[ ) provided 
several sets of EOS tables of nuclea r matter including hyperons. In this study, we utilize the 
EOS tables by llshizuka et al.l (120081 ) to examine both of attractive and repulsive cases for the 
undetermined S-iV inte raction supplementing the repulsive case reported in the Letter article 
( ISumiyoshi et al.ll2009l ). These Ishizuka EOS tables are based on an SUf{3) ext ended rela 



tivistic mean field (RMF) model and constructed as an extension of the EOS by lShen et al. 



( ll998al Jbl). which is widely utilized in various astrophysical simulations so far. Moreover, 
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Ishizuka et al.l (|2008[ ) prepared the hyperon EOS tables including the contributions of ther- 
mal pions for the extreme case where an effective mass of pions is equal to their rest mass 
in vacuum. We newly investigate the cases where both of hyperons and pions appear in 
failed supernovae, whil e the emergence of pions in nucleonic matter was already discussed 
jNakazato et aliboiObh . 



Therefore, the purpose of this paper is to report the detailed behavior of the black hole 
formation and the neutrino emission with the EOS with hyperons and/or pions. We as- 
sess especially the dependences on the hyperon-nucleon interaction and pion appearance in 
hyperonic matter for the black hole formation due to failed supernovae. We perform numer- 
ical simulations of the gravitational collapse, the c ore bounce, and the followin g evolution of 
proto-neutron star up to the black hole formation (jSumiyoshi et al.ll2007l . 120081 ). The general 
relativistic neutrino radiation hydrodynamics code, which solves the Boltzmann equations for 
neutrinos together with the Lagrangian hydrodynamics u nder spherical symmetry, is utilized 



to compute the dynarn ics as well as the neutrino signals (lYamadalll997l: lYamada et al 



1999 



Sumiyoshi et al.ll2005l ). The progenitor model with 4OM0 by IWoosley fc Weaven (119951 ) is 
adopted as the initial condition for the dynamical simulations. 

This paper is organized as follows. In § [2] we briefly describe our core collapse model 
including the numerics and initial setting. In §|3]we explain the hyperonic EOS. The results 
are discussed in § |U where we compare the dynamical features and neutrino signals with 
different hyperonic EOS's. §E]is devoted to summary and discussion. 



Numerical Simulations 



The numerical simulations are performed with the general relativistic neutrino radiation 
hydrodynamics code, which solves the Boltzmann equat ions for neutrinos together with the 



Lagrangian hydrodyn amics under spherical symmetry (lYamadal 119971 : lYamada et al.l Il999 



Sumiyoshi et al.ll2005l ). We consider four species of neutrino, z/g, ^^J. and t'^, assuming 



that the distribution function of {ut) is equal to that of z/^ 



A detailed description of 



the numerical simulations such as general relativistic hydrodynamics, transport and reaction 
rates of neutrinos, and resolutions can be found in ISumiyoshi et al.l (120071 ). We follow the 



dynamics from the onset of gravitational collapse of a progenitor through the core bounce 
and the post-bounce evolution of proto-neutron star by the accretion of outer layer, up to 
the formation of black hole. Note that we identify the b lack hole formation by finding the 
apparent horizon, as explained in iNakazato et al.l (120061 ). We evaluate the neutrino fluxes 
and spectra in detail up to the black hole formation owing to the exact treatment of neutrino 
transfer with a fully implicit method together with hydrodynamics. Note that, the neutrino- 
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hyperon reactions are omitted for simplicity, which will be revisited later. 



As an initial model, we employ the progenitor model of a AOMq star of lWoosley &: Weaver 

(Il995l ). This model contains an iron core of 2.0Mq. We use the profile of its central part up 



to S.OMq. While massive stars are suggested to lose their mass during the quasistatic evo- 
lutions, mass loss is not taken into account in the progenitor model we adopted. The mass 
loss rate is still uncertain in the theory of stellar evolution. The mass loss may affect dynam- 
i cs and neutrino signa l s of failed supernov ae through the density profile of the outer layer 



(Sumivoshi et al.ll2008l: iFischer et al. 



20091 ) as well as convection. Recently, lO'Connor fc Ott 



(I2OIII ) performed the core collapse simulations for several sets of the progenitor model with 
simplified neutrino transfer to find that the outcome depends severely on the mass loss pre- 
scription. Since our goal in this study is to investigate the influence of hyperons on the black 
hole formation, we fix the initial condition to a single progenitor model. 



3. Equation of State for Hyperon Matter 



The main input for the numerical simulations of the core collapse of massive stars is 
the EOS of dense matter in addition to the neutrino reaction rates because the dynamical 
time scale is long enough for particles other than neutrinos to equilibrate. One of the most 
widely used EOS tables is based on a Skyrme type (non-relativistic) mean field and the 
liquid-drop m odel of nuclei (jLattimer fc Swestylll99ll ). Another EOS table is based on an 



RMF model (ISugahara fc Toki 



1994), and nuclear formation is described in the Thomas- 



Fermi approximation ( IShen et al.lll998al Jbl). Hyperons are not included in these EOS tables, 
so that their applicable range may be limited to relatively low temperature and density 
regions, where hyperons do not emerge abundantly. Recently, several sets of EOS tables 
of nuclear matter including hyperons using an SUf{3) extended RMF model with a wide 
range of density, temperature, and ch arge fraction for nuin erical simulations of high energy 
astrophysical processes are presented (llshizuka et al.ll2008l ). 



We adopt the Ishizuka EOS in this study. This set of EOS is constructed as an extension 
of Shen EOS and smoothly connected with Shen EOS at low densities. In order to include 
hyperons, they take into account the potential depths of hyperons in symmetric nuclear 
matter at saturation density po, which are suggested in recent hypernuclear experiments. The 
potential depth of the A is well estimated as Uj^\pQ) ~ —30 MeV from the single-particle 
energies of many A hypernuclei. From the recently observed quasi-free S production spectra, 
it is considered that S hyperons would feel a repuls ive potential in nuclear matter, (po 
+30 MeV (e.g., 
potential value 



iarada fc Hirabavashill2005l . l2006l ). Nevertheless, having uncertainties of the 



Ishizuka et al. J2008h provided several sets of EOS tables using U, 



(TV) 



(Po 
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—30, +30, +90 MeV. Also for S hyperons, the an alyses of the twin hype rnuclear formation 
( lAoki et al.l Il995l ) and the S production spectra (IKhaustov et al.l |2000[ ) favor a potential 

depth of around U^^\po) 15 MeV. Note that more accurate measurement for the S 

potential is expected in the experiment at J-PARC. 

The S and S hyperons are particularly important in neutron stars, since nuclear matter 
can take a large energy gain from neutron Fermi energy and symmetry energy by replacing, 
for example, two neutrons with a proton and a negatively charged hyperon (S~ or 
If we adopt attractive potential for E hyperons, S~ would be the first hyperon to appear 
in neutron stars. With the recently favored potential around +30 MeV, E hyperons feel 
a repulsive interaction at high densities, and, hence, cannot appear until the baryon mass 
density pb 

around ~ 7 x 10^"^ g cm~'^ (see Figure 1 of llshizuka et al 



2 X 10^^ g cm ^. Instead of S hy perons, S hyperons are found to appear at 
7 X 10^"^ g cm~^ 



20081 ) . Therefore, we expect 



that the choice of hyperon interactions have an impact also on the dynamics and neutrino 
signal of failed supernovae. 

In the sets of Ishizuka EOS, the EOS tables including thermal pions are also constructed 
together with hyperon mixture. The thermal pions are treated in the minimum model where 
their effective mass is assumed to be equal to their rest mass in vacuum. This approach 
includes a simple s-wave Bose-Einstein condensation of pions , which is different f rom the 
pion condensation derived from p-wave vrA^ interaction (e.g., iKunihiro et al.lll993l ). Since 



s-wave 7iN interaction is repulsive, their effective mass becomes larger than that in vacuum. 
In this case, the pion population is suppressed. Thus, this EOS set corresponds to an 
extreme case where pions are overproduced, provided that p-wave nN attraction is omitted 
(lOhnishi et al.l 120091 ). We adopt this set to assess the maximum effect from thermal pions 
in hyperon mixture. 

Utilizing the set of Ishizuka EOS described above, we explore the influence of hyperon 
paying attention to the hyperon nucleon interaction and pion mixture. We investigate both 
of attractive and repulsive cases for the S potential. We perform the numerical simulation 
adopting the EOS sets with the potential depths (J/a, f/s, f/n) = (-30 MeV, +30 MeV, 
—15 MeV) for the repulsive case and (—30 MeV, —30 MeV, —15 MeV) for the attractive 
case. We perform also the corresponding cases adopting the hyperon EOS with pions. In the 
following, we refer to the repulsive EOS with pions, repulsive EOS without pions, attractive 
EOS with pions and attractive EOS without pions as RP, R, AP and A, r espectively. Note 



that results for the model R were already reported as the Letter article (ISumiyoshi et al. 



20091 ). For comparison, we also show the results for the EOS with pions and without hyperons 
and Shen EOS (purely nucleonic model), which are referred as NP and N, respectively. The 
model N was studied by ISumiyoshi et al.l (120071 ) in detail and the model NP was shown also 
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in iNakazato et al.l ( l2010bl ). where pions are treated with the same method as in Ishizaka 
EOS. 

The maximum mass of neutron stars gets generally lower due to new hadronic degrees 
of freedom. Hyperons are no exception. In fact, while the maximum masses of models N 
and NP are 2.2Mq and 2.OM0, respectively, those of every hyperonic models (RP, R, AP 
and A) are around I.QMq. Recen tly, the mass of the b inary millisecond pulsar J1614-2230 
was evaluated as 1.97 ± O.O4M0 (jPemorest et al.ll2010l ). This remarkable precision thanks 
to a strong Shapiro delay signature excludes the hyperonic models we adopted. While 
many calculations of the zero-temperature EOS for nuclear matter including hyperons were 
performed in various approaches, almost all of them cannot satisfy such a large value of the 
maximum mass. At present, Ishizuka EOS is a very limited model including hyperons which 
is available for astrophysical numerical simulations. We utilize Ishizuka EOS in this paper, 
where our focus is on investigating systematic differences due to the hyperon interaction 
and pion emergence. We believe that the systematics discussed in this study holds for more 
sophisticated EOS's which will be hopefully free from the maximum mass problem. Note 
that, ver y recently, Shen EOS was updated by the authors and A hyperon is taken into 
account (jShen et al.l |20 111 ). Their EOS does not include S and S, and its maximum mass is 



4. Numerical Results 

To begin with, we describe an outline of the core collapse evolution toward the black hole 
formation using the hyperon EOS. The radial trajectories of mass elements in the model with 
EOS R are shown in Figure [H The collapse is followed by the core bounce due to the nuclear 
repulsion as in ordinary core collapse supernovae. We measure the time from the timing of 
the core bounce thereafter. The central density at the bounce (3.2 x 10^^ g cm~^) does not 
differ among the models very much because the contribution of hyperons and pions is not 
dominant at this time as revisited later. The shock wave is launched due to the bounce. The 
propagation of the shock wave is also illustrated in Figure [H We can see that the shock wave 
stalls within 100 msec and recedes toward the surface of the central object, proto-neutron 
star under intense accretion. In the meantime, the proto-neutron star contracts gradually 
due to the increase of mass. It recollapses to the black hole abruptly at the critical mass. 
These pictures are the same for all the computed models. 

In Figure [21 we compare the time profiles of the central baryon mass density for all 
six models. The bounce corresponds to the spikes at t = and the black hole formation 
corresponds to the blow-ups. We can recognize that the time interval between the bounce 
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and black hole formation gets shorter as we put additional degrees of freedom, hyperons and 
pions. The time interval is 682 msec, 653 msec, 587 msec, 575 msec, 1345 msec and 1145 msec 
for the models with EOS's R, RP, A, AP, N and NP, respectively. The black hole formation 
occurs when the mass of proto-neutron star exceeds its critical mass due to the accretion 
of outer layer. Although the proto-neutron star is hot and lepton-rich, this critical mass 
is related with the maximum mass of neutron stars. As already mentioned, new hadronic 
degrees of freedom decrease the maximum mass of neutron stars. On the other hand, the 
mass accretion rate does not differ among the EOS models because hyperons and pions do 
not appear in the outer layer with low density and temperature. Hence the inclusion of 
hyperons and pions simply hasten the black hole formation. In contrast, the initial location 
of the apparent horizon, which is typically I.I-I.2M0 in the baryon mass coordinate, does 
not depend on the EOS. 

The neutrino duration time during black hole formation of the model with EOS R is 
~15% longer than that of the model with EOS A. Since the maximum masses of EOS's R 
and A do not differ very much, the duration time is more sensitive to the hyperonic matter 
EOS than the neutron star maximum mass is. Here we examine the profiles of the models 
with EOS's R and A to discuss the differences due to the hyperon interaction in the black 
hole formation. The profiles of key quantities at the selected times are shown in Figure [31 
At the core bounce, we can see that the hyperon potential does not make any difference 
because hyperons do not emerge yet. This is also evident from Figure HI where the profiles 
of particle fractions are shown. After the bounce, the density and temperature rise due to 
the contraction. This process is adiabatic as can be recognized from the entropy profiles. 
The electron-type lepton fraction does not also vary very much during this stage. 

Comparing the profiles at 500 msec after bounce, we can see that density and temper- 
ature of the model with EOS A are somewhat higher than that of EOS R while locations 
of the shock are the same for both models (see the velocity profile). As seen in Figure HJ A 
and S~ hyperons are populated near the center for the model with EOS A. By contrast, for 
the model with EOS R, the appearance of S~ hyperons is suppressed due to the repulsive 
potential. Therefore the density and temperature get higher for the attractive case because 
the EOS is more softened by the effect of hyperons. This is also the reason why the time 
interval between the bounce and black hole formation is shorter for the model with EOS A 
comparing with the model with EOS R. 

We discuss here the profiles at 2 msec before the black hole formation. Note that this 
moment corresponds to 680 msec after bounce for EOS R and 585 msec after bounce for 
EOS A. The central temperature of EOS A is lower than that of EOS R whereas their central 
densities are almost the same. This feature can be also recognized from Figure O where the 
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evolutions of the central densities and temperatures are shown. As already mentioned, 
evolutionary tracks are roughly isentropic lines. Therefore EOS R has higher temperature 
than EOS A for the same value of the entropy. This is equivalent to the fact that EOS A 
has higher entropy than EOS R for the same value of the temperature. In the attractive 
case, since E hyperons appear more easily and the number of particle species increases, the 
entropy gets higher for the fixed temperature. 

The peak of temperature profile resides not at the center but at the medium region 
(~O.7M0 of the baryon mass coordinate) for both models. This is because the entropy 
becomes high due to the shock heating at the medium region. The hyperon fractions are 
also different between the central and medium regions. For the repulsive case, S~ is the 
dominant negatively charged hyperon at the center, as is the case for cold neutron stars. 
In the medium region, the fraction of is larger than that of S~ because of the thermal 
population. On the other hand, for the attractive case, S~ appears more than S~ in the 
entire region. The abundance of negatively charged hyperons is also reflected in the charge 
chemical potential, which is the difference between the chemical potential of neutrons and 
that of protons, fin — fip- For the attractive case, the charge chemical potential is reduced 
owing to the emergence of E~. 

We now turn to the effect of pions. As shown in Figure [H qualitative features of the 
key quantities for EOS's RP and AP are similar to their counterparts without pions. Since 
the EOS is softened by the contribution of pions, the compression is accelerated both of the 
repulsive and attractive cases. As a consequence, the difference between the models with 
EOS's RP and AP is more minor than that of EOS's R and A. The appearance of pions results 
in the decrease of the electron fraction because the negative charge is shared by electrons 
and vr^. An influence of pions is more clearly seen in the profiles of the charge chemical 
potential. The chemical potential of charged pions, /Itt-, is equivalent to the charge chemical 
potential due to the balance, n ^ p + 7t~ . On the other hand, since the pions are thermally 
populated in our model, /Itt- is limited by their rest mass in vacuum, m^r- = 135 MeV. We 
can see //tt- saturates near 0.2Mq from the profile of the charge chemical potential of the 
model with EOS RP at 2 msec before black hole formation. 

The profiles of particle fractions with pions are shown in Figure [71 While there are 
quantitative differences between the models with and without pions, the pion mixing does 
not change the order of fractions for each hyperon. The fraction of vr^ at the center is 
larger for 500 ms after bounce than that for 2 msec before black hole formation because 
pions are replaced by the negatively charged hyperons for high density at 2 msec. During 
the contraction, pions appear so as to decrease the charge chemical potential. On the other 
hand, the octet baryons should be populated equally in high density or temperature limit. 
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where the charge neutrahty is aheady satisfied and pions are not necessary. Therefore, the 
existence of pions is hmited within the intermediate density regime. 

The core collapse of massive stars is accompanied by the emission of a large amount of 
neutrinos. The time profiles of luminosities and average energies of emitted neutrinos are 
shown in Figure |8] for all models. Note that and have the same type of reactions, the 
difference in coupling constants is minor and Ur (z^t) is assumed to be the same as [i^fj,). 
Therefore, we collectively denote these four species as z/^. The average energy presented here 
is defined by the root mean square value. The luminosity is corrected by the gravitational red- 
shift, though it is minor up to the end of our simulation at the apparent horizon formation. 
At the moment of black hole formation, the neutrino luminosity decreas es steeply due to 



the red-shift within a sho rt period and the emission stops suddenly (e.g.. iBaumgarte et al. 



19961 : ISeacom et allboOlH ^. Therefore the time mterval between the bounce and black hole 



formation almost corresponds to the duration of neutrino emission. As seen in the figure, 
the duration of neutrino emission is different among the models. This fact implies that the 
neutrino signal as an observable reflects the difference in the hyperon-nucleon interaction 
through the black hole formation. 

The neutrino signal of the stellar collapse can be seen as follows. The neutrino optical 
depth d{r) at the radius r is deflned as 

where Rg is the stellar radius and /mfp is the mean free path of neutrino, and the radius of 
the neutrino sphere Ri, is deflned as 

d{Ru) = I (2) 

for neutrinos with a typical energy. Roughly speaking, neutrinos are emitted from the 
neutrino sphere. In our cases, hyperons and pions appear only for the central region where 
the density and temperature are high enough for neutrinos to be trapped, i.e., d{r) ^ 1. In 
fact, the time evolutions of the neutrino signals are quite similar among the models except 
for the end points. We can recognize that the neutrino luminosity and average energy are 
slightly higher for the models with hyperons and/or pions. This is because they affect the 
location and temperature of the neutrino sphere through the softening of EOS. 

The neutrino luminosity summed over all species is given approximately by the accretion 



^However see also ISekiguchi fc Shibatal (|201ir ) for moderately rotating collapse 
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luminosity (iThompson et al.ll2003l ) 



B ' ^'^) 

where G, M and are the gravitational constant, the mass accretion rate and the mass 
enclosed by R^, respectively. In our case, as is the case for ordinary supernovae, the proto- 
neutron star is formed before the black hole formation, and neutrinos are emitted mainly 
on the surface of the proto-neutron star. Thus we can regard Ry and as the radius and 
mass of the proto-neutron star, respectively. The radius of the proto-neutron star becomes 
smaller for the models with hyperons and/or pions because of the softening of EOS. Thus 
the neutrino luminosity gets higher. On the other hand, the average energy of neutrinos 
is approximately proportional to the temperature of the neutrino sphere, T^. When R^ 
decreases because of the softening of EOS, increases. This is the reason why the average 
energy of neutrinos gets higher for the models with hyperons and/or pions. These features 
are remarkable for the case with an attractive E-A^ interaction because E~ hyperons are easy 
to appear and help the softening of EOS. Incidentally, in order to examine these differences 
in the neutrino signals using actual events of failed supernoyae, mor e detailed statistical 



analyses is necessary based on the method in iNakazato et al.l (l2010al ). Progenitors in our 
Galaxy will be favorable for the statistical studies. 

We have checked that the omission of neutrino-hyperon reactions in our simulation is a 
reasonable approximation through the evaluation of diffusion time scale. In fact, hyperons 
appear only in the central region and neutrinos in the central region take time to diffuse out. 
Thus the absence of the neutrino-hyperon reactions will not make a difference for the early 
phase. The diffusion time from the radius r can be evaluated as 

rRs g^/^^/ 

^<'-* = / S7^' 

where c is the velocity of light. In Figure HI we compare the diffusion time of z/g with 30 MeV, 
which is typical energy of emitted neutrinos, of the model R for the ca ses with and withou t 



the neutrino-hyperon reactions based on the reaction rates shown in iReddy et al.l (1l998l ). 
We can recognize that the difference is subtle for 500 msec after bounce because the hyperon 
fraction is tiny. Since the difference resides inside IM©, where the diffusion time is >500 ms, 
the influence of hyperons will be seen in the neutrino signals after >500 ms disregarding the 
evolution of the mean free path. At 2 msec before black hole formation, neutrino-hyperon 
reactions makes a noticeable difference inside 1.3Mq, but the diffusion time is still long with 
r > 250 ms. Therefore we can conclude that the omission of the neutrino-hyperon reactions 
does not affect the results on the emitted neutrino signals. 
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5. Summary and Discussion 

In this paper, we have presented the emergence of hyperons in black-hole-forming failed 
supernovae with the focus on the hyperon-nucleon interaction and the pion mixture. A series 
of core-collapse simulation for a progenitor model with AOMq have been performed solving 
general relativistic neutrino radiation hydrodynamics under spherical symmetry. We have 
adopted a few sets of EOS tables of nuclear matter including hyperons using an SUf{3) 
extended RMF model. In particular, we have paid attention to the difference in the S 
potential, whether attractive or repulsive. The appearance of thermal pions have been also 
taken into account including a simple s-wave Bose-Einstein condensation. The collapse leads 
to the core bounce and the birth of proto-neutron stars with increasing mass, while the shock 
wave stalls and the explosion fails being different from ordinary core-collapse supernovae. 
We have found that the inclusion of hyperons and pions promotes the recoUapse towards 
the black hole formation by terminating the proto-neutron star epoch because the EOS is 
softened due to new hadronic degrees of freedom. Since S hyperons appear more easily, the 
softening is more efficient for the attractive S potential cases comparing with the repulsive 
counterparts. We have found also that thermal pions act as an additional source of softening 
together with hyperons. 

We have also evaluated the amount of neutrinos emitted from the failed supernovae. 
For the spherically symmetric nonrotating collapse, the time interval from the bounce to the 
black hole formation corresponds to the duration of neutrino emission. The softer EOS's 
with new hadronic degrees of freedom lead to the smaller critical masses and the earlier 
timing of black hole formation. Therefore the duration of neutrino signal is shorter for the 
models with hyperons and/or pions. The time evolutions of the neutrino signals before the 
black hole formation are found rather similar among the models other than the difference in 
duration. This is because they appear only inside the proto-neutron star, which is opaque 
to neutrinos. However, there are slight differences in the rise of energies and luminosities 
around the end points in the duration of neutrino signal. They would provide a clue to assess 
the physics of strangeness in addition to the information from the time duration. In fact, 
the event number of neutrinos at SuperKamiokande detector is evaluated for the G alactic 



events as ~10 , which could be a target of statistical studies (iNakazato et al.ll2010al ) 



As already mentioned, there are large uncertainties in the mass loss prescription of the 
progenitor model. It affects the s tellar structure and, therefore, the time interval between the 



bounce and black hole formation ( Sumiyoshi et al.ll2008l : lFischer et al.ll2009l : IO'Connor fc Ott 



2OIII ). Thus it is our concern that the difference in the progenitor mass produces changes 
in the duration of the neutrino signal as the difference in the EOS does. However, there 
may be the observational chances for probing the progenitor apart from the neutrino. For 



- 13 - 



instance, since A ve can determine th e direction of the progenitor to some extent by the neu- 
trino detection f lAndo fc Satdl2002l ). optical follow-up observations may be possible. While 
the neutrino sphere is swallowed by the black hole for in ~1 sec, the free fall time at the 
stellar surface (photosphere) is 1 / \/Gp ~ 10^ sec. By observing the stellar disappearance, 
we may be able to get more information on the stellar structure. Combining such additional 
data, we believe that the assessment of the EOS from the neutrino signal is promising. 

As a final note, we would like to emphasize that the knowledge about the properties 
of dense matter at extreme conditions is essential for the studies on the astrophysical black 
hole formation such as not only the collapse of massive stars but also the merger of binary 
neutron stars (ISekiguchi et al.ll201ll ) and so on. 
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0.2 0.4 0.6 0.8 

time (s) 

Fig. 1. — Radial trajectories of mass elements as a function of time after bounce in the model 
with EOS R. The trajectories are plotted for each O.O2M0 in baryon mass coordinates. Thick 
solid lines denote the trajectories for 0.5, 1.0, 1.5, 2.0 and 2.5Mq. The location of the shock 
wave is shown by a thick dashed line. 
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Fig. 2. — Time profiles of tlie central baryon mass density. Thick solid, thin solid, thick 
dashed, thin dashed, thick dot-dashed and thin dot-dashed lines correspond to the results 
for EOS's R, RP, A, AP, N and NP, respectively. 
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Fig. 3. — Profiles of the density (upper left), temperature (upper center), electron-type 
lepton fraction (upper right), velocity (lower left), entropy per baryon (lower center) and 
charge chemical potential /i„ — /Xp (lower right) at bounce, 500 msec after bounce and 2 msec 
before black hole formation. The results for EOS's R and A are shown by thin solid and thick 
dashed lines, respectively. Note that the 2 msec before black hole formation corresponds to 
680 msec after bounce for EOS R and 585 msec after bounce for EOS A. 
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Fig. 4. — Profiles of the particle fractions at bounce (left), 500 msec after bounce (center) 
and 2 msec before black hole formation (right). The upper and lower panels show results for 
EOS's R and A, respectively. 



-20- 



50 



40 



> 



30 - 



0) 



10 












5x10^4 

density (g/ cm^) 



10 



15 



Fig. 5. — Evolutions of the central density and temperature (solid lines) with hyperon 
fraction contours for the electron fraction Ye = 0.3 (dashed lines). Thick and thin lines 
correspond to EOS's R and A, respectively. The dot-dashed line shows the trajectory of the 
model with EOS N. 




Fig. 6. — Same as Figure |3] but for EOS's RP (thin solid lines) and AP (thick dashed lines). 
Note that the 2 msec before black hole formation corresponds to 651 msec after bounce for 
EOS RP and 573 msec after bounce for EOS AP. 
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Fig. 7.- 



Same as Figure H] but for EOS's RP (upper panels) and AP (lower panels). 
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Fig. 8. — Luminosities (upper plots) and average energies (lower plots) of the emitted neu- 
trinos as a function of time after bounce. The panels correspond, from left to right, to v^) 
z/g and (= ^'^t, '^t, ^^i^ ^r)- The results for the models with EOS's R (thick dashed lines), 
A (thick solid lines) and N (thin solid lines) are shown in the top panels, and the results for 
the models with EOS's RP (thick dashed lines), AP (thick solid lines) and NP (thin solid 
lines) are shown in the bottom panels. Vertical lines represent the end point of the neutrino 
emission. 
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Fig. 9. — Profiles of tlie diffusion time of witli 30 MeV for tlie model R at 500 msec after 
bounce (dashed lines) and 680 msec after bounce (solid lines). Thick lines represent cases 
with the neutrino-hyperon reactions while thin lines represent cases without the neutrino- 
hyperon reactions. 



